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ABSTRACT 

Be stars make up almost 20% of the B star population, and are rapidly rotating stars surrounded by a disc; however the origin 
of this rotation remains unclear. Mass transfer within close binaries provides the leading hypothesis, with previous detections 
of stripped companions to Be stars supporting this. Here, we exploit the exquisite astrometric precision of Gaia to carry out the 
largest to date comparative study into the binarity of matched samples of nearby B and Be stars from the Bright Star Catalogue. 
By utilising new “proper motion anomaly” values, derived from Gaia DR2 and DR3 astrometric data alongside previous values 
calculated using Hipparcos and Gaia data, and the Gaia provided RUWE, we demonstrate that we can identify unresolved binaries 
down to separations of 0.02”. Using these measures, we find that the binary fractions of B and Be stars are similar between 0.04 
— 10”, but the Be binary fraction is significantly lower than that of the B stars for separations below 0.04”. As the separation 
range of these “missing” binaries is too large for mass transfer, and stripped companions are not retrieved by these measures, 
we suggest the companions migrate inwards via binary hardening within a triple system. This confirms statistically for the first 
time the hypothesis that binary interaction causes the Be phenomenon, with migration causing the dearth of Be binaries between 
0.02 - 0.04’’. Furthermore, we suggest that triplicity plays a vital role in this migration, and thus in the formation of Be stars as 
a whole. 
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1 INTRODUCTION cloud, the rotation thus a result of inherited angular momentum 
(Bodenheimer 1995); the envelope of the star could be spun up via 
momentum transport from a contracting core (Granada et al. 2013; 
Hastings et al. 2020); or, crucially for this study, it could be a re- 
sult of mass and angular momentum transfer via binary interaction 
(Shao & Li 2014; Bodensteiner et al. 2020a). 


If binary interaction is the source of this rapid rotation, then one 


Be stars make up around 20% of the total population of B-type stars 
(Bodensteiner et al. 2020a) and are defined, in part, by the presence 
of Balmer emission lines. These lines are thought to originate from a 
slowly outflowing ionised circumstellar gas-disc (Porter & Rivinius 
2003; Rivinius et al. 2013; Jones et al. 2022). While the dynamics of 
these discs have been studied extensively (Wisniewski et al. 2010; ; pas oe 
Draper et al. 2014; Suffak et al. 2022), the mechanism by which could expect that the binary statistics of B and Be stars will differ. 
they form is undetermined. Consensus suggests that the formation This line of reasoning provided the basis for the study performed m 
is related to the rapid rotation observed in practically the entire Oudmaijer & Parr (2010), which found the binary fraction of B and 


Be population (Zorec et al. 2016, 2017), and which is suggested to Be-type stars to be similar between separations of 0.1 and 8” based 
be at or near the critical velocities of such stars (Townsend et al. on a sample of 36 B and 37 Be stars, and concluded that binarity, 


2004; Granada et al. 2013). Such rapid rotation could thus allow for at least at these scales, could not be responsible for all Be stars. 
processes such as turbulence (Townsend et al. 2004) or non-radial However, while a number of studies into the binarity of Be stars exist 
pulsations (Owocki & Cranmer 2002; Semaan et al. 2018) to con- (Boubert & Evans 2018; Klement et al. 2019; Hastings et al. 2021 
tribute significantly to the formation of the decretion discs around and El-Badry et al. 2022, among others) and suggest the stripping of 
Be stars due to decreased effective gravity at the surface of such stars the companion to the Be star as the cause (Han et al. 2002; Nazé et al. 


(Granada et al. 2013). This, however, only leads to a further question 2022; El-Badry & Quataert 2021; El Badry et al. 2022), very few 
- what exactly causes Be stars to rapidly rotate? comparative studies into the binarity of both B and Be-type stars 


Three mechanisms have been suggested as the origin of this char- exist. The next most recent points of comparison are Abt & Levy 
acteristic rotation; it could be a relic from the parent molecular (1978) and Abt & Cardona (1984), each from around 40 years ago, 
and which similarly conclude that binarity cannot be responsible for 


the Be phenomenon at the scales they probe. However, Abt & Levy 
* E-mail: py17jd@leeds.ac.uk (1978), which utilises radial velocity variations to detect binarity, 
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suffers from the effects of small sample sizes, while Abt & Cardona 
(1984) is hampered by inhomogeneous biases as it is, in effect, a 
review of the then available literature. 

With the advent of high-precision astrometric all-sky surveys, 
we are now able to mitigate these issues and investigate the mul- 
tiplicity of large samples. Recent studies have employed astrometric 
data from surveys such as Hipparcos (Van Leeuwen 2007) and Gaia 
(Gaia Collaboration et al. 2021) in order to detect binarity in the 
overall stellar population, both in the case of wide (El-Badry et al. 
2021) and close (Kervella et al. 2019, 2022) binaries. The time is 
thus ripe for such an investigation into the binarity of B and Be stars. 

Here we utilise the Proper Motion Anomaly (PMa), as introduced 
in Kervella et al. (2019) (see also Brandt 2018), which is a mea- 
sure of the difference between long and short term proper motions 
of a point source, in order to determine the fractions of B and Be 
type stars in visually unresolved binary systems. Alongside previ- 
ously determined PMa values from Kervella et al. (2022), calculated 
using Hipparcos and Gaia DR3 data (a similar catalogue can also 
be found in Brandt 2021), we use new values, calculated for the 
first time in this work using a combination of Gaia DR2 and DR3 
data (Gaia Collaboration et al. 2018, 2022). Using these metrics we 
probe binary systems with separations between 0.02 and 1.1”, and 
with magnitude differences up to 4mag. 

Additionally we utilise the Gaia-provided Renormalised Unit 
Weight Error (RUWE) to further constrain the binary fractions of 
these two populations (Stassun & Torres 2021). Thus we are able 
to bridge the gap in binary separations between, and expand upon 
previous radial velocity (Abt & Levy 1978) and adaptive optics 
(Oudmaijer & Parr 2010) studies of these stars. Additionally, by 
querying a sample over an order of magnitude larger than present 
in the most recent similar study (Oudmaijer & Parr 2010), this paper 
provides the largest scale comparative study of B and Be star binarity 
to date. 

This paper is structured as follows: in Section 2 we detail the 
two binary detection methods used (the PMa and the RUWE), and 
demonstrate their sensitivities and limitations; in Section 3 we report 
the results of utilising these methods on samples of B and Be stars; 
and in Section 4 we discuss the implications of these newly reported 
binary fractions on models for the formation of Be stars. 


2 METHODS 
2.1 Binary Detections via Proper Motion Anomaly 


For any binary system we can define two similar variables, the centre 
of mass, which the two component stars orbit, and the photocentre, 
which describes the location of the apparent point source in an un- 
resolved system. In general, the positions of these two points differ, 
thus the photocentre also orbits the centre of mass - leading to a 
deviation from the linear proper motion expected for single sources 
in the case of unresolved binary systems. 

This deviation from the expected linear proper motion, referred to 
as the Proper Motion Anomaly (PMa), Ay, can be detected through 
a comparison of long and short term proper motion measurements 
(Uit and ust respectively), with the long term effectively tracing the 
motion of the centre of mass and the short term tracing the sum of 
the centre of mass and relative photocentre motions. An illustration 
of this principle can be found in Figure 1. The PMa can thus be 
expressed as: 


AH = Hst — Hit (1) 
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Figure 1. Cartoon illustrating the principle behind the Proper Motion 
Anomaly. The large and small stars are the primary and secondary com- 
panion respectively, the black dot is the centre of mass and the yellow dot 
the photocentre. Two sets of observations are taken at A and B, with a time 
interval of At between them. ua and upg are the short term proper motions 
measured during their respective observations, and zag is the long term 
proper motion, determined from the change of observed position over At. Ar- 
rows labelled upc are the instantaneous proper motions of the photocentre 
in the frame which the centre of mass is stationary. 


Thanks to current advances in the precision of proper motion 
measurements, alongside high precision positional measurements 
allowing for long term proper motions to be calculated, the PMa 
thus provides a powerful tool in the search for unresolved binary 
candidates within large samples of stars. 

This measure is sensitive to the period of a given binary system. 
Periods significantly longer than the time elapsed over the long term 
proper motion measurement, At will have a deviation in proper mo- 
tion too small to be detected. Similarly, systems wherein a whole 
number of complete orbits have occurred during the measurement 
of the long term proper motion will not be detected, again as the 
deviation in long and short term proper motions will be too small to 
detect. 

Following this logic, peak sensitivity of the PMa should be 
achieved for systems where an i orbital periods have elapsed be- 
tween the two observations, as the absolute change in orbital phase 
is at a maximum here. However, in reality, due to the time-window 
smearing effects present over the observation periods of the two con- 
stituent astrometric surveys (Kervella et al. 2022) there is a general 
decrease in sensitivity towards smaller periods, with peak sensitiv- 
ity being achieved for periods of P ~ At (Kervella et al. 2022). This 
also explains why, although of much higher astrometric precision, the 
DR2-DR3 derived PMa’s probe larger periods than the Hipparcos- 
DR3 PMa values (as will be shown in the following sub-subsection). 
Below we will probe the sensitivity limits of the various methods 
both theoretically and empirically. 


2.1.1 Separation Limits of the Proper Motion Anomaly 


Whereas the sensitivity of the PMa to mass as a function of orbital 
period (and thus linear separation) has previously been investigated 
(Kervella et al. 2022), the sensitivity to angular separation remains 
unexamined. As detailed in Appendix A, using an idealised theoreti- 
cal model, consisting of a face-on, circularly orbiting binary system, 
we arrive at a minimum detectable separation of 0.01” for PMa 
values calculated using Hipparcos and Gaia DR3 data. Additionally, 


when considering the PMa calculated using Gaia DR2 and DR3 data, 
the theoretical minimum detectable separation is 0.0025”. However, 
this computation does not take into account the effect of the time- 
window smearing during the PM measurement - which leads to the 
DR2-DR3 PMa losing out on the shorter separations. 

It is thus useful to obtain an empirical value as a point of compar- 
ison. In order to evaluate this sensitivity of the PMa to varying sepa- 
rations, we use a sample of binary systems listed in the Washington 
Double Star (WDS, Mason et al. 2022) catalogue. The WDS contains 
acompilation of multiple star systems for which at least one measured 
separation is available within the literature. The sources in the WDS 
were cross-matched with the list provided in Kervella et al. (2022) 
using the xMatch service provided by CDS, Strasbourg (Pineau et al. 
2020). No restrictions on mass or spectral type were placed on the 
stars within the binaries in order to generate the largest, general sam- 
ple. Reported higher order multiple systems were then removed for 
the sake of simplicity, as were any systems with separations greater 
than 10”. While the PMa primarily probes the unresolved regime of 
separations by nature, this 10” limit was selected in order to ensure 
that no unpredictable effects were materialising at higher separations. 

For the resulting 11 116 unique binary systems, we obtained values 
of the PMa from the Kervella et al. (2022) catalogue (which uses 
data reported in the Hipparcos re-reduction (Van Leeuwen 2007) 
and the Gaia DR3 catalogues, giving Ar ~ 25 yrs) and calculated 
new values, using data available in the Gaia Data Release 2 and 
Data Release 3 (Gaia Collaboration et al. 2018, 2021; At ~ 0.5 yrs). 
The methodology used in Kervella et al. (2022) was utilised for the 
calculation of these new values, with calibration performed via re- 
calculating Hipparcos-Gaia DR3 PMa’s replicating the reported PMa 
values exactly and reported errors to within a precision of 1%. 

The final sample contains 10 283 sources with reported PMa val- 
ues; note that not all values surpass the binary detection threshold of 
S/N>3. The majority of these systems had separations reported to a 
precision of 0.1”. 316 systems have separations reported as 0” in the 
WDS Catalogue (i.e. their separations were below 0.1%). For each 
of these systems, higher precision separations were retrieved from 
supplemental WDS material (Mason et al. 2022). Where orbital so- 
lutions had previously been reported, the separation available at the 
Gaia DR3 epoch was used. 

Figure 2 shows the median values of the signal-to-noise ratio 
for the Hipparcos-Gaia DR3 PMa for the WDS sample at varying 
separations. These systems have been placed into bins of width 0.1” 
at separations greater than or equal to 0.1’’. At separations below 
0.1” the systems have been placed into 10 bins such that each bin 
has an equal number of systems within it (approximately 30 per bin). 
For these bins the mean reported separation of the systems within it 
is plotted. 

Based on the data presented in Figure 2, it can be seen that for 
the DR3-Hipparcos PMa, there is an increase of SNR towards lower 
reported separations. On average, this PMa remains above the detec- 
tion threshold of 3 for all separations between 0.02 and 0.7” (based 
on the median values). As a result, 0.02 — 0.7” appears to be an 
empirical range of separations detectable via this PMa; at least half 
of all WDS binary systems with separations within this range are 
detected (note that we wouldn’t expect for all binaries within this 
range to be detected due to projection effects). This is in agreement 
with the calculated theoretical minimum detectable separation and 
the expectation that the upper limit occurs somewhere around the 
spatial resolution limit of the surveys. 

It is also to be noted that high PMa systems at separations greater 
than the resolvable limit (i.e. with separations greater than ~ 1’’) exist 
in non-negligible numbers. These systems have reported separations 
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Figure 2. Median values of the Hipparcos - Gaia DR3 and Gaia DR2 - 
DR3 PMa S/N’s and the RUWE, of systems reported as binary in the WDS 
(Mason et al. 2022). Here values greater than 3 constitute a detection of 
binarity by either PMa, and values greater than 1.4 constitute a detection by 
the RUWE - see Section 2.2 for details. 


much greater than the PMa can detect (as it measures the motion of 
an essentially unresolved binary), and it is thus likely that they are, 
in fact, as of yet undetected higher order multiple systems. 

Alongside the examination of the Kervella et al. (2022) reported 
Hipparcos-Gaia DR3 PMa values, an equivalent consideration of 
Gaia DR2-DR3 PMa values has been conducted. As there does not 
exist an equivalent catalogue as in Kervella et al. (2022), these PMa’s 
were thus calculated via the same process, but using data from the 
Gaia Data Releases 2 and 3. 

The median values of these PMa SNR values at a given separation 
can also be found in Figure 2. Again, a rise in the median PMa can 
be seen towards lower separations, here beginning at separations of 
approximately 1.1’’, increasing towards a peak at around 0.6”, and 
remaining above the detection threshold down to 0.2”. While this 
lower limit is in disagreement with the theoretical limit of 0.0025’ 
derived in Appendix A, the theoretical considerations largely ignore 
the effect of smearing over the observing timeframe. As the DR2- 
DR3 PMa is smeared across the entirety of the long-term proper 
motions timeframe, this will mean that the lowest separations are 
not probed, although the number of potential observable objects 
increases significantly (as Gaia has observed orders or magnitude 
more objects than Hipparcos). Additionally the increase in precision 
achieved since Hipparcos has both resulted in the upper limit of 
separations detectable to be expanded beyond the Gaia resolution 
limit, and increased the total fraction of systems detected within the 
sensitive limit. 

These results remain consistent when examining B and Be type 
binaries within the WDS (Mason et al. 2022), with systems within 
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Figure 3. Median PMa of WDS systems, as a function of the difference in 
magnitude of the two stars for the Hipparcos-Gaia DR3 and Gaia DR2-DR3 
PMa’s and the RUWE. Both plots have been cropped to a maximum difference 
in magnitude of 8mag, as the number of systems with magnitude differences 
greater than this decreases to the point where meaningful conclusions cannot 
be drawn. Detections are reported out to differences of up to 12mag in all 
three cases. 


the two sensitive ranges being detected at significantly higher rates 
than those outside that range. 

We thus show that the Hipparcos-Gaia DR3 PMa is sensitive to 
binary systems with angular separations between 0.02—0.7’’, and the 
Gaia DR2 - DR3 PMa is sensitive between separations of 0.2 — 1.1”. 


2.1.2 Magnitude Limits of the PMa 


In order to evaluate the limitations of the Proper Motion Anomaly, 
we must also consider the sensitivity of this measure to the difference 
in magnitude between the two components. Here systems contained 
within the WDS catalogue are once again considered, now limiting 
the sample to those systems with reported separations between 0 and 
1” in line with the more favourably sensitive range determined in the 
previous section. 

In the case of the Hipparcos-Gaia DR3 PMa, and as can be seen 
in Figure 3, examining the difference in the WDS reported apparent 
magnitude reveals that the PMa shows some preferential sensitivity 
within the 1.0 to 3.5mag range, peaking at around AVinag = 3mag. 

In the case of the Gaia DR2-DR3 PMa, an equivalent examina- 
tion reveals a general increase in sensitivity across most magnitude 
differences available via the WDS sample, with the preferentially 
detected region expanding to systems with a difference in magnitude 
0 < AVinag < 4.5, peaking at around AVmag = 2.5, although sig- 
nificant numbers of detections are still made at greater differences in 
magnitude. 
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2.2 Binary Detections via the RUWE 


The Renormalised Unit Weight Error, or the RUWE, is a reduced 
chi-squared parameter that describes the ‘goodness-of-fit’ of a single 
body astrometric model to an observed source. As such, it would 
follow that such a parameter is sensitive to the presence of unre- 
solved binary companions - with the photocentre wobble caused by 
such a companion causing deviations from what is expected of a 
single body, resulting in a worse fit astrometric solution. Indeed, 
the RUWE parameter has been suggested to act as a binarity in- 
dicator (Belokurov et al. 2020), with values significantly above the 
normalised value of 1 (i.e > 1.4 as per Stassun & Torres 2021) being 
regarded as indicative of the source being ‘badly fit’ by the astromet- 
ric solution and thus potentially binary. 

Like the PMa, the RUWE varies in sensitivity according to the 
period of a binary system, with the RUWE having some function of 
sensitivity to binarity determined both by the period of the system 
in question and by the time baseline over which the RUWE is deter- 
mined (approximately six months in the case of Gaia DR2 and DR3), 
with Penoyre et al. (2022) suggesting that for Gaia DR3, RUWE is 
an effective measure of binarity for periods of around a month to a 
decade. 

However, due to sensitivity to other (non-multiplicity related) fac- 
tors which can inflate the RUWE (Belokurov et al. 2020), it is impor- 
tant to note that there may be some fraction of false positive detections 
- for example if the sample contains stars known to have circumstellar 
discs or extended regions of emission (Oudmaijer et al. 2022). This 
is unlikely to hamper our study of B and Be stars however, as the 
discs of Be stars are very small when compared to the resolution 
achieved by Gaia DR3 (Quirrenbach et al. 1997; Wheelwright et al. 
2012). 


2.2.1 Separation Limits of the RUWE 


The separation limits of the RUWE were determined using the same 
process as for the PMa above. Here, the detection threshold is given 
to be 1.4 (as suggested by Stassun & Torres 2021). Based on the data 
presented in Figure 2, we find that the RUWE exhibits a sharp and 
marked increase above the detection threshold at binary separations 
below the resolution limit of Gaia (0.7’’). This increase peaks at 
approximately 0.09”, and, on average, remains above the threshold 
down to separations of 0.04”. Thus, we find that the RUWE is 
sensitive to those binary systems with separations ranging between 
0.04” and 0.7”. 


2.2.2 Magnitude Limits of the RUWE 


As with the PMa, it is useful to quantify the sensitivity of the RUWE 
to the magnitude differences of the binary systems probed. Figure 
3 demonstrates, by repeating the same process as used to determine 
the magnitude differences the PMa is sensitive to, that it would seem 
that the RUWE is primarily sensitive to systems with magnitude 
differences between 1.5 and 3.0mag, peaking somewhere between 
2.0 and 2.5mag. 

We have thus shown that, via these astrometric methods, we can 
consistently detect visually unresolved binaries with separations as 
low as 0.04” using the RUWE, and 0.02” using the Hipparcos-Gaia 
PMa and with magnitude differences up to 4.5mag, depending on 
the measure used. Thus, thanks to the precision achieved by modern 
astrometric surveys, we will be able to probe the origins of Be stars 
with a statistical basis heretofore unavailable. 


Table 1. Populations of main sequence B and Be type stars within the sample, 
divided by spectral type. 


Spectral Type B Be 


BOV 12 4 
BIV 47 10 
B2V 104 31 
B3V 97 19 
B4v 31 12 
BSV 111 13 
B6V 43 8 
B7V 48 4 
B8V 156 13 
B9V 258 9 
Total 907 | 123 


3 THE BINARITY OF B AND BE STARS 


3.1 Sample Selection 


Here we utilise the same initial sample of B and Be stars ex- 
tracted from the Bright Star Catalogue (BSC) 5th Revised Edi- 
tion (Hoffleit & Warren 1995; Hoffleit & Saladyga 1997) as used in 
Radley et al. (in preparation). This sample was then cross-matched 
within 5” against the Gaia Data Release 3 (Gaia Collaboration et al. 
2021) data. Stars without proper motion measurements in the Hippar- 
cos or Gaia Data Releases were removed. Stars with non-V luminosity 
classes were also removed in order to ensure that the comparisons 
made between our two samples are between like sources that are sub- 
ject to the same selection effects. A summary of the spectral types of 
the stars included in this sample can be found in Table 1. 

This sample provides an almost thirty-fold increase in the number 
of B-type stars, and a four-fold increase in the number of Be-type 
stars compared to similar comparative studies (Abt & Levy 1978; 
Oudmaijer & Parr 2010) that probe a well defined separation range. 
Additionally, by virtue of being selected from the same catalogue 
of the brightest stars in the sky, these stars will be subject to the 
same selection effects, and have similar brightnesses and parallaxes, 
thus ensuring their usefulness in comparing how B and Be type stars 
differ in their astrometric properties. These stars are also amongst the 
most nearby ones, thus allowing us to detect the physically closest 
companions whilst remaining within the range of angular separations 
detectable by our methods (see Sections 2.1.1 and 2.2.1), with the 
lower end of our detectable range, 0.02” equating to 5 au at the mean 
distance to the stars in these samples. 


3.2 Detection of Unresolved Binaries 


We first consider the complete, unmodified samples of 907 B and 
123 Be stars, utilising a combination of PMa values found within 
the catalogue contained within Kervella et al. (2022), calculated us- 
ing Hipparcos and Gaia DR3 data, and new PMa values calculated 
using data reported in Gaia Data Releases 2 and 3, as well as the 
Gaia-provided RUWE. These measures allow us to detect unresolved 
binary systems via differences in long and short term proper motions 
(for the PMa; see Section 2.1) and via deviations from a single body 
astrometric fit (in the case of the RUWE; see Section 2.2). 

Using the Hipparcos-Gaia DR3 PMa, which probes the range 0.02 
— 0.7” (see Section 2.1.1), equating to around 5 — 180 au at the 
average distance of 260pc, we arrive at binary fractions of 42 + 2% 
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and 28 + 4% for the B and Be samples respectively!. This proves to 
be a greater than 30 difference. 

In contrast, the Gaia DR2 - DR3 PMa values reveal similar B and 
Be binary fractions in the 0.2 — 1.1” range (around 50 — 290 au at 
260pc). In this case, the fraction of B stars detected as being in a 
binary system is found to be 27 + 1% and the equivalent for Be stars 
is found to be 29 + 4%. 

Finally, the fractions of B and Be stars detected via applying a de- 
tection threshold of RUWE > 1.4, tracing binaries with separations 
between 0.04 — 1.1” (10 — 290 au), were afound to be 27 + 1% and 
20 + 4% respectively. See Table 2 further down for a summary of 
these binary fractions. 


3.3 Correcting Biases 


Although the global properties of the B and Be stars in the BSC 
are very similar, the two samples exhibit slightly different properties 
across a number of observables. For example, by virtue of being 
rarer, the Be stars are on average further away than the B stars, with 
the average B star in our sample sitting at a parallax of 5.26mas with 
a standard deviation of 3.38mas, whereas the average Be star sits at 
a parallax of 3.81mas with a standard deviation of 2.42mas. Like- 
wise the average apparent magnitudes of the B and Be populations 
are 6.06mag (standard deviation of 1.07mag) and 5.95mag (standard 
deviation of 1.32mag) respectively. In addition to this, the two pop- 
ulations differ in the distribution of their spectral types (for example, 
almost 30% of the B stars are B9V, while only 7% of the Be stars are 
B9Ve). As such, the variations in the distributions of these properties 
must be taken into account to ensure that the binary fractions are not 
affected by them. 

Thus, in order to remove the effects of these biases, a variety of 
sub-samples of B stars were created. These included random samples, 
alongside “‘semi-random” samples wherein the B population was 
binned according to a given variable, mentioned above, and stars were 
randomly picked from each bin in order to match the distribution of 
said variable in the Be sample. In each case a Kolmogorov-Smirnov 
test was performed in order to confirm that the sub-sample was likely 
drawn from the same distribution as the Be population. The binary 
fraction of each sub-sample was then determined using the PMa and 
RUWE detection thresholds detailed above. 

A summary of these samples and their binary fractions can be 
found in Figure 4. Therein, the greater than 30 dearth of Be binaries, 
when compared to the number of B binaries, is still seen to be present 
at the tightest separations probed by our methods. 


3.4 Wider Separation Binaries 


While the focus of this study is on close binaries, we now consider 
the presence of binary companions outside the 0.02 — 0.7” range. 
Oudmaijer & Parr (2010), who probe separations of 0.1 — 8”, deter- 
mined binary fractions of B and Be stars at 29 + 8 and 30 + 8%. This 
finding could imply that the difference in the two binary fractions 
disappears towards larger separations. 

Exploiting the added information provided by Gaia, we performed 
a basic search for spatially resolved B and Be binaries by cross- 
matching the sample depicted in Fig. 4, and looking for sources 
within a radius of 10’’. Any ‘duplicate’ sources that were retrieved 


' Uncertainties here are calculated using the binomial method. Error bars are 
given at lo. 
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Figure 4. Summary of the binary fractions of B and Be stars, as determined by Hipparcos-Gaia DR3, Gaia DR2 - DR3 PMa values, Gaia DR3 RUWE values 


and resolved binaries (see text). Also presented are the binary fractions of the 


various sub-samples of the B type stars, taken in order to account for any biasing 


effects. ‘Matched’ here indicates that a sub-sample of the B stars has been created such that the distribution of a given variable in the sub-sample is drawn from 
the same parent distribution as the Be population. Note that the ‘Wide Binaries’ plot is noisier due to smaller number statistics. The Hipparcos - Gaia DR3 


values, probing the smallest separations, show a significantly lower Be binary 


frequency (in each case a greater than 30 difference), while all other approaches 


return similar binary fractions to within the errors. Based upon this we report a dearth of Be binaries with separations between 0.02 and 0.7”, and, as discussed 
in Section 3.5 we further constrain these missing binaries to between 0.02 and 0.2”. 


were then deemed to be probable companions to the primary if they 
satisfied: 


|Wprimary = duplicate | < Oprimary,w + Tduplicate,w (2) 
where w is the parallax of either star, and Cw is the uncertainty 
corresponding to either parallax. 

Based on this basic test, 6+ 1% of B and 10+3% of Be V-type stars 
are found to have potential resolved companions within a distance 
of 10”. These values remained consistent even when randomly sam- 
pling or when generating semi-random samples of B stars intended 
to match the spectral type, parallax or apparent magnitude distribu- 
tions of the Be stars. A summary of these results can also be found in 
Figure 4. In the case of each of the distribution matching samples, a 
Kolmogorov-Smirnoy test was performed in order to verify that the 
distribution of each given parameter in the B star sample is likely 
being drawn from the same parent distribution as the Be star sample. 

Presuming that each of these additional sources is indeed a binary 
companion, poses an interesting question, that being why do Be stars 
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have binary companions at a lower rate than their B type equivalents 
at close (i.e. < 0.7” separations, but at a similar rates at greater 
separations? 


3.5 Further Constraints on Binary Separation 


In this study we have determined, through the use of various astro- 
metric measures, between separations of 0.02 — 0.7”, the binary 
fraction of Be stars is significantly lower than that of their B star 
equivalents, while they are similar at larger separations. Let us con- 
tinue by further identifying at which separation the B binary stars 
become more prevalent. 

Table 2 reveals that by comparing the two PMa values and looking 
at those systems detected by the Hipparcos-Gaia DR3 PMa but not by 
the Gaia DR2-DR3 PMa, we can constrain the lack of Be binaries to 
between 0.02 and 0.2”. We also note that at separations greater than 
0.2”, B and Be stars have similar fractions of their populations in 
binary systems - in agreement with the previous most recent compar- 


Difference in B and Be star binarity 7 


Table 2. Percentage of B and Be stars detected as binary at varying ranges of separations. By utilising combinations of the previously explored astrometric 
detection methods we are able to constrain the greater than 30 death of Be binaries to between 0.02 and 0.2”. 


Detection Method 


Separation Range /” B Be 
0.02 < x < 0.2 28 + 1% 17+ 3% 
0.02 < x < 0.7 42+2% | 28+4% 
0.2<x<1.1 27+1% | 29+4% 
0.04 < x < 0.7 27+1% | 20+4% 

O.1<x <8 29 + 8% | 30+ 8% 
0.7< x < 10 6 + 1% 10 + 3% 


ative study of these stars (Oudmaijer & Parr 2010) Additionally, 
we note that the RUWE detects statistically similar rates of B and 
Be stars as being in binary systems. Thus, by virtue of the RUWE’s 
sensitivity stretching between 0.04 and 0.7’’, we are able to confine 
the dearth of Be binaries even further - restricting them to between 
0.02 and 0.04”. 


4 DISCUSSION 


In previous sections we have established, for the first time, ranges of 
angular separation and magnitude difference within which the PMa 
and RUWE are sensitive to binarity. 

We have then used these measures to reveal a greater than 30 
dearth of Be stars in binary systems at separations unresolved by Gaia, 
compared to B stars, we then go to further constrain this absence of 
Be binaries to separations between 0.02 and 0.04”. 

These results pose a conundrum, as one might naively assume that 
this implies that the Be phenomenon more favourably occurs via sin- 
gle star pathways, rather than the currently favoured binary pathways. 
Therefore, in order to continue to pursue the binary pathway for Be 
formation, we must consider: what sort of companion will remain 
undetected by these methods, whilst also remaining consistent with 
such models of Be formation? 

In this section we will posit some explanations for this difference 
in binary fraction, alongside the implication for models regarding the 
formation of the Be phenomenon. 


4.1 Stripped Companions 


One such solution invokes the presence of a companion whose enve- 
lope has been stripped by the putative Be star. Accretion from such a 
companion would provide a source for the initial spin-up mechanism 
that eventually results in the formation of the Be stars’ characteristic 
disc (El-Badry et al. 2022; Jones et al. 2022), while also allowing the 
companion to become low-mass, dim or close enough to the host star 
(Frost et al. 2022) to be unlikely to be detected via the PMa. 


2 Although we do not wish to make quantitative comparisons between values 
reported by our different detection methods, a larger binary fraction detected 
by the Hipparcos-Gaia DR3 PMa (42 + 2%) than the Gaia DR2-DR3 PMa 
(27+ 1%) would be expected in the case of a log-flat distribution of separations 
(as per Opiks law; Opik 1924). This is because the separation range probed 
by the Hipparcos-Gaia DR3 PMa (0.02 to 0.7”) is almost twice the width 
in log-space as that of the Gaia DR2-DR3 PMa (0.2 to 1.1”). Conversely, 
by themselves the values reported for the Be stars may indicate a distribution 
decreasing towards to closer separations that we are sensitive to. 


Detected by the Hipparcos-DR3 
PMa but not the DR2-DR3 PMa. 
Detected by the Hipparcos-DR3 
PMa. 

Detected by the Gaia DR2-DR3 
PMa. 

Detected by the RUWE. 


Detected by Oudmaijer & Parr (2010) 


Resolved as binary by Gaia. 


Such a model has been suggested to be the route through which 
20 — 100% of field Be stars form (El-Badry & Quataert 2021), 
and a significant number of these stripped companions to Be 
stars have been reported in the literature (El-Badry et al. 2022; 
Bodensteiner et al. 2020b; Frost et al. 2022; see Table 3). Addition- 
ally, no close main sequence companions to Be stars have been 
confirmed (Bodensteiner et al. 2020a), which is consistent with what 
would be expected in the case of a stripped companion origin for 
Be stars (although it should be noted that this study only considers 
massive Be stars - i.e. spectral types earlier than B1.5). As such, they 
may be responsible for hydrogen-poor supernovae and neutron star 
mergers resulting in gravitational wave sources. 

The periods proposed to be required in order for stripping to begin 
by Han et al. (2002) are of order 0.1 — 10 days for Common Envelope 
Ejection and 400 — 1500 days for Roche Lobe Overflow. Using the 
average distance to the Be stars in our sample (~ 260pc), and as- 
suming a circular, Keplerian orbit, these stars would require a binary 
companion at a separation of 0.02” (approximately 5 au) or closer 
in order for stripping to initiate. While, the process of mass transfer 
would then act to widen the orbit of this companion, the resultant 
stripped star would be low-mass and faint. We would thus expect that 
the majority of stripped companions, if they indeed exist, would be 
missed by both our PMa values and the RUWE as they will too close 
to the host Be star, or too low mass and dim to be detected. 

In support of this, PMa and RUWE values were obtained for 16 Be 
stars confirmed via X-ray observations to host a stripped companion 
(all of those from Nazé et al. 2022). In order for the above suggestion 
to have any validity, these stars, that are known to host stripped 
companions, should not be detected as binaries using the measures 
detailed within this paper. A summary of these objects and their Gaia 
data can be found in Table 3. 

Seven of these 16 systems are detected by either method; 5 by the 
Hipparcos - Gaia DR3 PMa, and 5 by the RUWE, with 3 objects de- 
tected by both. However, 5 of these 7 detected systems (the exceptions 
being HD 194335 and HD 200310) have known additional compan- 
ions (Mason et al. 2022) at separations within the range the PMa’s 
and RUWE are sensitive to. In addition, these companions are at 
separations significantly larger than required for the periods reported 
by Nazé et al. (2022), and are therefore inconsistent with being the 
reported stripped companion. An undetected, third, companion has 
also been suggested in the case of HD 200310 to explain the eccentric 
orbit of the stripped companion (Klement et al. 2022), making the 
system a potential triple. Thus, the true fraction of detected stripped 
companions would seem to be either 6 or 12% depending on how 
one wishes to treat the case of HD 200310. Therefore, we reach the 
conclusion that the stripped companions listed in Nazé et al. (2022) 
are by and large not detected by the PMa or RUWE, whose detec- 
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Table 3. Summary of Gaia DR3 - Hipparcos PMa and RUWE values for Be 
stars quoted in Nazé et al. (2022) as being known to host stripped compan- 
ions. PMa and RUWE values marked with an asterisk (*) constitute binary 
detections. 


Name Spectral Type | PMa 


RUWE 


g Per B1.5Ve 0.23 | 1.804 
HD 29441 | B2.5Vne 1.17 1.125 
HD 41335 | B1.5IV-Ve 8.13* | 1.515* 
HD 43544 | B3Ve 0.66 1.047 
HD 51354 | B3Ve 0.79 1.055 
HD 58978 | B0.5IVe 2.52 | 0.941 
HD 60855 | B2Ve 1.52 1.097 
HD 113120 | B2IVne 6.93* | 2.385% 

x! Aps B2Vnpe 1.48 1.060 
HD 152478 | B3Vnpe 1.82 | 0.990 
HD 157832 | B1.5Ve 1.60 | 0.735 

28 Cyg B3IVe 1.55 1.113 
HD 194335 | B2IlIle 3.81* | 1.041 

59 Cyg BlVe 6.38* | 2.724* 

60 Cyg B1Ve 6.31* | 0.142 

8 Lac A BlIVe 0.87 | 1.480* 


tions turn out to trace a third body in these systems instead. Thus, 
the potential for the Be phenomenon being caused by the presence 
of stripped companions remains consistent with our results. 


4.2 Triples 


Intriguingly, the separation range of the “missing” Be binaries (0.02 
- 0.04”; with 0.02” corresponding to Sau) is typically too large for 
mass transfer to have occurred (typically less than 4au is required per 
the periods reported in Han et al. 2002). The natural next question 
to ask is then, why would the Be stars lack binaries in this particular 
separation range, while the stripping of a close companion should 
occur at smaller separations? 

A clue can be found in the fact that the systems listed in Table 3 
show a high incidence rate of higher order multiplicity, with 56+ 18% 
of this sample of Be stars known to host a stripped companion having 
additional companions (Mason et al. 2022) at separations larger than 
those reported in Nazé et al. (2022). Additionally this number is also 
somewhat elevated with the expected incidence rate of higher or- 
der multiplicity posited by Moe & Di Stefano (2017) - being around 
40%. 

The presence of such a third companion can neatly explain the 
absence of detected Be binaries at larger separations, while also al- 
lowing for any companion stars to become close enough for mass 
transfer to occur; it is well known that higher order multiplicity can 
result in the hardening of an inner binary. Indeed, a third body in- 
creases the occurrences of migration and eventual binary interactions 
significantly (Toonen et al. 2022; Kummer et al. 2023). 

Recent models (Preece et al. 2022) suggest that of the systems 
that eventually form stripped companions - half form via binary 
pathways and half form via triple pathways. This can take place by 
either mass transfer occurring within the inner binary - with the outer 
tertiary star potentially becoming unbound - or a merger of the inner 
binary and mass transfer occurring between the two remaining stars 
(Preece et al. 2022). 

Our results therefore imply that close binary interactions are re- 
sponsible for the formation of Be stars, which constitute 20% of the 
B star population, with a migratory effect causing the lack of Be 
binaries in the 0.02 to 0.04”. Moreover, we suggest that triplicity 
must play a vital role in triggering this migration, and thus in the 
formation of Be stars as a whole. 
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4.3 Potential Limitations of the Methods 


We will now address the various potential limitations of the methods 
as described in Sections 2 and 3. 


4.3.1 Sample Selection 


First, we note that the Be phenomenon is known to be transient, 
with the disc known to disappear and reappear on short timescales, 
which can lead to the potential misclassification of Be stars as B stars 
Rivinius et al. (2013). Thus, contamination of the B star sample with 
misclassified Be stars could effect the previously detected binary 
fractions for the B star sample. However, as the Be sample should 
have low contamination (as observation of the emission line features 
associated with the Be phenomenon is required for classification as a 
Be star, even if these features later disappear), the detected Be binary 
fractions should be representative of the Be population as a whole 
- and thus representative of any misclassified Be stars. Therefore, 
should there be a significant amount of Be contaminants in the B star 
sample, then they will act to bring the true B binary fraction down - 
thus the true B binary fraction would differ more significantly from 
the Be fraction. 

We also note that, as this sample is selected from 
the BSC (Hoffleit & Jaschek 1982; Hoffleit & Warren 1995; 
Hoffleit & Saladyga 1997), these stars are necessarily bright. Whilst 
Gaia astrometric fits become systematic dominated at G < 13 
(Lindegren 2020), the resultant increase in uncertainties associated 
with this effect ensures that the PMa is still a viable metric - as a 
signal-to-noise ratio is used to quantify binary detection. Addition- 
ally, in the calculation of the new Gaia DR2 - DR3 PMa’s correc- 
tions to the respective proper motions of bright stars as proscribed 
by Lindegren et al. (2018) and Cantat-Gaudin & Brandt (2021) have 
been applied. Steps have also been taken to ensure that the two sam- 
ples are subject to the same selection effects, and these are detailed in 
Section 3.3, thus any systematic effects should be consistent between 
the two samples. To check this, all stars brighter than Smag were 
removed from the B and Be samples; which resulted in zero change 
to the reported dearth of Be binaries within the 0.02 — 0.2” range. 


4.3.2 Variability 


While photometric variability has not been quantitatively considered 
here, we anticipate that it is unlikely to effect the PMa unless the 
variability is occurring both asymmetrically and on a scale similar 
to the sensitive separation range, here we look at this in more detail. 

First, consider a single star exhibiting symmetrical variation; in 
this case the position of the photocentre does not deviate from the 
centre of mass - thus leading to no proper motion anomaly. When 
we then introduce an unresolved companion, the position of the 
photocentre therefore deviates from the position of the centre of mass, 
however the magnitude of the deviation varies with the variability 
cycle. Thus, while such variability shouldn’t produce a false positive 
detection of binarity, it is possible for it to induce a false negative. 
The likelihood of this false negative thus depend on two factors: the 
magnitude of the variability, and the timescale on which it occurs 
(as observations taken at the same point in the variability cycle will 
obviously be unaffected, and variability with a period shorter than 
either individual observing period will be smeared, diminishing the 
effect on the PMa). 

In the case of asymmetric variability - i.e. that caused by phe- 
nomena such as mass ejection, or variations in surface brightness in 
stars such as red super-giants (Chiavassa et al. 2022), or asymmetry 


within a disc (Meilland et al. 2007; Stee et al. 2013) - deviation of 
the photocentre from the centre of mass will occur. However, such 
asymmetry would have to be on a similar scale to the PMa’s sensitive 
range to in order for such deviation to be construed as a false positive 
detection of binarity. 

In the specific case of this study, we anticipate little effect of the 
photometric variability of Be stars (Labadie-Bartz et al. 2017, 2022) 
on our PMa values. This is due to the observed periods of Be vari- 
ability being significantly shorter than the observing periods of Gaia 
DR2, DR3 and Hipparcos (668, 1038 and 1227 days respectively), 
with the longest period variability found by Labadie-Bartz et al. 
(2017) being approximately 200 days. Similarly, variability caused 
by mass loss onto the Be disc is unlikely to significantly effect the 
PMa, as the discs of Be stars, with sizes of order milli-arcseconds, are 
very small when compared to the resolution achieved by Gaia DR3 
(i.e. 0.7’; Wheelwright et al. (2012); Quirrenbach et al. (1997)). 
Additionally, by similar logic, asymmetry of the discs themselves 
(Meilland et al. 2007; Stee et al. 2013) is also expected to have little 
effect on the PMa measurements of these stars. Likewise, the stars 
are necessariliy smaller than their disc, they are expected to be too 
small to have asymmetric surface brightness that would effect the 
PMa. 


5 CONCLUSIONS 


In this paper we have performed the largest scale comparative study 
of B and Be-star binarity to date, utilising both the Hipparcos - Gaia 
DR3 PMa values found in Kervella et al. (2022) and new Gaia DR2 
- DR3 PMa values calculated for the first time here, alongside the 
Gaia provided RUWE. 

We evaluate the limits in which these methods are sensitive to 
binarity. Particularly we determine for the first time the range of 
angular separations they are capable of reliably detecting, through 
the use of WDS (Mason et al. 2022) data. We thus find that the 
Hipparcos - Gaia DR3 and Gaia DR2 - DR3 PMa’s are sensitive 
to binaries with separations in the 0.02 — 0.7” and 0.2 — 1.1” 
ranges respectively. Likewise the RUWE is found to be sensitive to 
separations between 0.04 and 0.7”. 

Through the use of these astrometric methods in combination 
with one another, we report a greater than 30 dearth of Be binaries 
compared to B stars at separations between 0.02 and 0.2”, with the 
fractions of B and Be stars found to be in binary systems at these 
separations being 28 + 1 and 17 + 3% respectively. 

We posit that this apparent lack of Be binaries provides evidence 
for the presence of stripped companions to the Be stars, with the 
mass-transfer process by which they are created having been pre- 
viously suggested as an origin for the Be phenomenon. Such stars 
would go undetected by the measures used within this paper, as they 
are too close, low-mass and dim to induce a significant enough de- 
viation of the photocentre from the centre of mass. We also find 
that Be stars previously reported to have a stripped companion have 
additional companions (i.e. triples and higher order multiples) at a 
rate somewhat elevated from the expected rate for B stars as a whole. 
Therefore, we suggest that migration of binary companions to Be 
stars to separations at which mass-transfer can occur, triggered by a 
third companion, must play an important role in the formation of Be 
stars. 
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APPENDIX A: DERIVATION OF THE THEORETICAL 
LIMITS OF THE PMA 


In order to obtain a theoretical, lower bound binary separation de- 
tectable via the PMa, it is easiest to consider a simple binary system. 
Take, for example, a system of two intermediate mass (i.e. 2 — 5Mo) 
main sequence stars in a circular orbit. In such a system, where both 
stars are of a known mass, the centre of mass is easily derivable. 
Likewise, by applying the mass-luminosity relationship: 


ales iy ee (Al) 


for such stars, as given by Griffiths et al. (1988), the location of 
the photocentre is similarly able to be determined. 

The minimum theoretically detectable separation will occur for a 
system where (a) the difference between the centre of mass and the 
photocentre is at a maximum and, (b) where the change in absolute 
orbital phase is at a maximum (i.e the observations are spaced half 
an orbit apart). 

The latter of these criteria (b) is fairly easy to evaluate, provided we 
know the masses of the two stars. However, we must first evaluate the 
former criteria (a) in order to determine the mass ratio at which the 
largest offset between the photocentre and centre of mass is achieved. 
This proves to be a somewhat more involved process. 

By normalising the separation between the two stars to 1, it is 
possible to write the photocentre and centre of mass equations as: 


41,M/L + 42,M/L = 1 (A2) 
M,a1,.m = M202,M (A3) 
Liai yp = Loa, (A4) 
ô = |a1,m -41,L\ (A5) 
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(2) companion from the centre of mass (M) or photocentre (L) de- 
pending on the subscript, M is the mass of an individual star, L is 
the corresponding luminosity and 6 is the offset between the centre 
of mass and the photocentre. 

From this set of equations, alongside the mass-luminosity relation 
given by Equation A1, it is possible to rewrite Equations A3 and 
A4 in terms of a1 m| /z and the mass ratio of the two stars, Mı /M2, 
giving the expression of 6 as: 


M (M2*° _ M2 (M2 
Mı Mı Mı Mı 


By taking the derivative of this with respect to the mass ratio, 
and solving for dô/d (%) = 0, it is possible to obtain the mass 
ratio for which the maximum offset between the centre of mass and 
photocentre is achieved. This is found to occur for a mass ratio, 
M/M; = 0.47 for intermediate mass main sequence stars. 

If we now assume a circular orbit with a stationary centre of mass, 
then the photocentre motion is resultant entirely from the orbital 
motion of the system. Thus, as the maximum change in proper motion 
occurs for observations separated by oe L orbital periods, where n 
is some integer greater or equal to zero, in this case resulting in 
the change in proper motion being twice the instantaneous proper 


é=| 1)7! +1)7!}] (A6) 


motion. Thus, we can simply express this change in proper motion, 
Apas: 


2nr 
Ap = 2st = wrô = p’ (A7) 


where ust is the short term proper motion measurement, which 
here can be reduced to the tangential velocity of the photocentre. P 
is the period of the orbit and r is the binary separation. 

While this approximation of the relation between the change in 
proper motion and the period of a binary system may seem to imply 
that the proper motion changes increasingly towards shorter period 
binaries, this approximation ignores the effects of smearing on the 
actual observationaly data. In reality this smearing across the ob- 
servational timeframe causes a drop in sensitivity at shorter periods 
(Kervella et al. 2022), with the maximum sensitivity being achieved 
for periods similar to twice that of the time elapsed between the 
two observations. Here, this resolves to a period of ~ 50yrs for the 
Hipparcos-Gaia DR3 PMa, and ~ lyr in the case of the Gaia DR2- 
DR3 PMa. 

Using this calculated PMa, in concert with standard errors for 
the three surveys in question (in this case the Hipparcos re-reduction 
(Van Leeuwen 2007) and the Gaia DR2 and 3 (Lindegren et al. 2018, 
2021) although other values can be used as and where appropriate), 
a lower limit of separation can be obtained. This is achieved via 
propagating the standard errors through the calculations required to 
determine the PMa (as detailed in Kervella et al. (2022)) and using 
this to generate a signal-to-noise ratio. As we have an expression for 
the change in proper motion that is a function of the separation we 
can thus write: 


“= 5 (A8) 


where N is the standard error of calculating the PMa. Thus by 
setting the above to equal the Kervella et al. (2022) defined detection 
threshold of 3, and solving for r we arrive at a theoretical minimum 
binary separation detectable by the PMa. In the case of the PMa 
obtained from Hipparcos re-reduction and Gaia DR3 data this lower 
limit is estimated to be 0.01”, and in the case of the Gaia DR2- 
DR3 PMa this is 0.0025”. However the effect of the observational 
smearing is underestimated, particularly in the case of the Gaia DR2- 
DR3 PMa, where the PM measurement period covers the entirety 
of long term proper motion’s timeframe - likely resulting in less 
sensitivity to small separations than predicted here. 


This paper has been typeset from a TEX/IATEX file prepared by the author. 
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